Hidden Dimer in the Frenkel-Kontorova Model 
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The incommensurability of the supercritical Frenkel-Kontorova model is decomposed into a family 
of dimer type "defects" by appropriate decimations. Interestingly, this hidden dimer results in 
Bloch-wave type excitations of the renormalized chain which appear in the disguised form of multi- 
step phonon modes for the original system. We call these intriguing excitations stepon modes and 
conjecture that they play a key role in determining the localization boundary in systems exhibiting 
Anderson localization. 
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PACS numbers: 05.45. +b,63. 50. +x,64.60.Ak 

In this paper, we describe a very interesting relation- 
ship between two seemingly different models that have 
been in the forefront of theoretical physics: the Frenkel- 
Kontorova (FK) model Q and the random dimer model 
H . This relationship establishes the existence of prop- 
agating Bloch-type phonons in the pinned phase of the 
FK model and resolves the mystery of novel multi-step 
excitations found recently [Q. 

The incommensurate FK model consists of a chain 
of balls connected by Hooke's springs in a periodic si- 
nusoidal potential of the strength K where the average 
spacing a between the balls is incommensurate with the 
periodicity of the potential. The equilibrium positions of 
the balls x n are given by the iterates of the standard map 



sin(27ra:„ ). 
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For small K, the iterates cir6 confined on an invariant 
circle of the standard map. The breakup of the invariant 
circle at the critical value K c corresponds to the pin- 
ning transition in the FK model accompanied with the 
disappearence of the zero frequency phonon mode. The 
phonon modes are linear excitations of the equilibrium 
configurations of the balls satisfying the equation 



(2) 



where the eigenvalue E is related to the phonon frequency 
u> as E — —lo 2 + 2, and e„ = K cos(27ra„) is the onsite 
phonon potential. 

At the so-called anti-integrable limit K — > 00 ||, the 
balls sit at the bottom of the potential wells and therefore 
e n is constant for all sites n. When K decreases, the limit- 
ing uniform phonon potential splits into more and more 
distinguishable levels. As shown in Fig. 1, this split- 
ting generates a symbolic representation of the lattice. It 
turns out that by expanding the potential upto any order 
1/K P the resulting incommensurate lattice with no ap- 
parent short range clustering can be decimated into a lat- 
tice with dimer and tetramer "defects" . As p increases, 



the decimated structure becomes more and more com- 
plex. The discovery of a hidden dimer helps in tracing the 
origin of intriguing "stepon" excitations with multi-step 
eigenf unctions that were found in our recent numerical 
study of the supercritical parameter region K > K c [Q. 
The central result of this paper is that the stepon modes 
of the FK model arc in fact propagating Bloch waves on 
the renormalized lattice. 

The models with dimer type defects attracted a great 
deal of attention a few years ago || due to the discov- 
ery of subtle derealization mechanisms which explained 
the conduction in polymers. The delocalized modes orig- 
inated from special resonance conditions where the re- 
flected waves from neighboring defect sites interfered de- 
structively resulting in the ballistic transport in the sys- 
tem. Detailed theoretical studies showed that a variety 
of extended defects with reflection symmetry could be 
made to exhibit this resonance by tuning a control pa- 
rameter so that the resonance condition was satisfied by 
an energy close to the Fermi energy. This turns out to 
be true also for the dimer and tetramer defects in the 
renormalized lattices of the supercritical FK model. 

For our analysis, we study a generalized version of the 
FK model with two independent parameters K and A so 
that in Eq. (|2|) e„ = Acos(27ra n ). The model reduces 
to the FK model when the two parameters are equal. 
This extension of parameter space is analogous to the 
analytic continuation into complex plane and will turn 
out to be crucial in understanding the stepon modes. 
In the 2-parameter model, the parameter A controls the 
strength while K determines the smoothness of the onsite 
"phonon" potential as the underlying invariant circle of 
the standard map undergoes the transition by breaking of 
analyticity at the critical value K — K c |J . In the region 
K < K c , the potential e„ is a smooth function of the ef- 
fective phase 9 = {ner + </>} where the brackets denote the 
fractional part. In this case the model falls into the uni- 
versality class of the famous Harper equation which 
exhibits a transition from extended to localized states. 
On the other hand, for K > K c the phonon potential is 
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determined by an underlying invariant Cantor set (can- 
torus) of the standard map. Therefore, the 2-parameter 
model provides one with a new class of quasiperiodic sys- 
tems that interpolates nicely between the Harper equa- 
tion [K = 0) and the FK phonon equation. As seen 
later, this important feature helps in understanding the 
absence of Anderson localization for the FK phonons. 

Fig. 2 shows the phase diagram of the model obtained 
using an exact decimation scheme discussed in our ear- 
lier papers The model exhibits an extended phase, 
a localized phase as well as a critical phase with self- 
similar eigenfunctions. A novel aspect of the model is 
the emergence of infinitely many curves from the "cor- 
ner" A, K — > oo along which the eigenfunctions are rep- 
resented by an infinite series of step functions of the ef- 
fective phase (see Fig. 3). We call these solutions 
stepon modes. All except the right most curve intersect 
the standard FK limit K — A. We referred to these in- 
tersections as degeneracy points in our earlier numerical 
work [0] because at these parameter values, the nontriv- 
ial scaling properties of critical phonons "degenerated" 
into the trivial scaling of the stepon modes. It is in- 
teresting to note that the degeneracy curves accumulate 
at the critical parameter value K = K c for the circle- 
cantorus transition of the standard map. The parameter 
region outside these curves consists of either the critical 
or localized phase. The localized phase for large values 
of A intertwines with the degeneracy curves. The fact 
that the localized modes reside inside the tongues of the 
degeneracy curves suggests that the stepon modes play 
a special role in the localization phenomenon in models 
where the underlying potential is not smooth. However, 
the specifics of that role remains eluded to us at present. 

The above phase diagram suggests a perturbative ap- 
proach to understand the origin of the stepon modes. We 
develop a perturbation theory near the anti-integrable 
limit by expanding the equilibrium configurations of the 
balls x n with k = 1/K as the expansion parameter. 
In the cantorus regime, the hull function X defined by 
x n = X(na + cj)) is a convergent series of step functions 
§. We write X{9) = X Q {6) + X x {6)k + X 2 {6)n 2 + 
where X o (0) = \ + Int{6). Substituting the above ex- 
pansion into Eq. ([!]), we generate a similar expansion for 
the "phonon" potential A cos[27rAf (#)]. Truncating this 
expansion at the order k p leads to a potential with a 
finite number of steps. This in turn leads to symbolic 
representation of the lattice of the type shown in Fig. 1, 
where the number of required symbols increases with the 
order of perturbation. 

For general p, the symbolic representation can be very 
complex. However, it turns out that by an appropri- 
ate decimation, the lattice can always be represented by 
one of the two possible forms which happen to be the 
ones found at p = 2 and p = 3. Different p cases are 
distinguished from each other by different renormalized 
onsite "energies" and coupling terms. In the following, 



we summarize the essential features of the symbolic rep- 
resentation. 

Each symbol (except the one for the potential min- 
imum) is associated with two symmetrically placed 6- 
intervals where the potential takes the corresponding 
constant value. For example, for p — 2 the symbol a is 
attached to the ^-intervals [0, c), [1— a, 1) and the symbol 
b to the interval [er, 1 — a). Asp increases, these intervals 
split into subintervals in a systematic way. The borders 
of intervals, which are discontinuities of the potential, are 
obtained as the first p — 1 forward and backward iterates 
of the map Oj+i = Oj + c (mod 1) with 8q = 0. The split- 
ting of the levels as p increases can be shown to follow two 
possible patterns depending upon whether the new dis- 
continuities land within the middle interval or not. If the 
two new discontinuities at the order p land inside order- 
(p— 1) intervals other than the middle one, both intervals 
are split into two parts whose relative lengths are deter- 
mined by the golden mean. We label the longer ones of 
these new subintervals by a and the shorter ones by b. 
The rest of the intervals can be labeled in an arbitrary 
way. This pattern emerges for p — 3, 4, 6, 7, .... The other 
possibility is that both new discontinuities land inside the 
middle interval which breaks up into three parts a, b, and 
a. In this case the total length of the a-subintervals is 
related to the length of the 6-subinterval by the square 
of the golden mean. This happens to be the case for 
p = 2,5,18,.... 

It turns out that the lattice sites which are not labeled 
as a or b can be decimated. In the renormalized lattice, 
the complexity of the original symbol dynamics mani- 
fests itself in the renormalization of the onsite energies 
as well as the coupling terms. Because of the symmetric 
alignment of the decimated blocks jl0| , the renormalized 
couplings between the a-type sites alternate between two 
values. Thus, the whole problem reduces to showing the 
existence of propagating wave solutions for the lattice 
with dimer and tetramer "defects" . 

It suffices to discuss the p = 3 perturbation the- 
ory in detail (p = 2 gives the trivial solution) as the 
other renormalized lattices have a similar form. The 
truncated potential in the p = 3 case takes three dis- 
tinct values e a , e&, and e c as shown in Fig. 1. Here 
e a = A(-l + 2tt 2 k 2 (1 - 6k)), e b = A(-l + 2tt 2 k 2 (1 - 8k)), 
and e c = —A. The decimation of the sites c results in 
the renormalized lattice ...aaaabbaabb... with the renor- 
malized onsite energy e a = (1 — {E — e a )(E — e c )), while 
the renormalized couplings between two neighboring a- 
sites alternate between V\ = 1 and V<z = (E — e c ). The 
decimated lattice can be shown to have a traveling wave 
solution provided the parameters satisfy the condition 

(e b - e a ){E - e a )(E - e c ) = (e b - e„) + (e c - e„). (3) 

Each decimated c-site gives rise to the phase shift ft given 
by 
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(e 6 - e a )e- lk + 1 



1 



(4) 



(eft - e a )e ik 

where E = £b + 2cos(k). Assuming that these phase shifts 
amount to the total phase shift f2(n) at the nth site of 
the renormalized lattice, the traveling wave solution can 
be expressed as e l l kn + n ( n )} . 

It should be noted that within the defects, the form of 
the wave is rather complicated. However, by decimating 
the sites labeled by c, the complexity has been absorbed 
into the renormalization of the onsite energies and the 
coupling terms. For k = 0, we always have 51 = except 
when e& — e a = —1 for which the k — solution corre- 
sponds to Q = 7T. This special solution is actually very 
similar to the one obtained by Dunlap et al. || in their 
study of the random dimer model. We would also like to 
point out that our general solution is valid for the type 
of defects we discuss irrespectively of whether they orig- 
inate from quasiperiodic, chaotic or correlated random 
processes. 

The existence of numerically obtained stepon modes 
of the type described in Fig. 3 implies that delocalized 
modes exist also for the system obtained by taking into 
account the full expansion upto infinite order. In this 
case, the stepon eigenfunction is represented by an infi- 
nite series of step functions. Therefore, the stepon modes 
of the FK model can be thought of as Bloch waves on an 
infinitely many times renormalized lattice. 

Another important result of our analysis is the fact 
that it predicts the asymptotic form of the degeneracy 
curves shown in Fig. 2. Using the explicit values of the 
onsite energies results in the following relationship be- 
tween the two parameters k and A for the three rightmost 
branches: 

— -- 3 + o(0 
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was obtained using the p — 3 perturbation theory 
t)) while the p = 4 theory gives all three solutions 
(|5j|7|). We conjecture that the rest of degeneracy curves 
can be explained by the higher order perturbation theory. 

In summary, our numerical results along with the dec- 
imation and the systematic perturbation theory demon- 
strate that the incommensurate FK model in the pinned 
phase is related to the class of models with corre- 
lated defects and exhibits propagating wave solutions on 
the renormalized lattice. Moreover, the extended two- 
parameter model is the first known case which exhibits 
both correlated defects and the localization transition. 
The propagating stepon modes seem to play a special 
role in determining the onset to Anderson localization. 
Proper understanding of how the stepon modes deter- 
mine the localization boundary may be difficult as the 



perturbation theory in any finite order does not pre- 
dict the existence of exponentially localized modes in the 
model. 

The study of the two-parameter model helps one to un- 
derstand an important distinction between the FK model 
and the Harper equation. First of all, the absence of lo- 
calization in the FK model is due to the fact that the 
nonlinear potential is not strong enough to localize the 
phonons. The extended FK model, where the strength 
and the smoothness are controlled independently, does 
exhibit the localization transition for A >> K. However, 
the localized phase for K > K c exhibits an important dis- 
tinction from that of the Harper equation. The localized 
eigenfunctions (once the exponentially decaying part is 
factorized out) for K < K c exhibit universal self-similar 
fluctuations characterized by a unique strong coupling 
renormalization fixed point |p[. In contrast, these fractal 
fluctuations in the cantorus regime K > K c appear irreg- 
ular and defy any simple renormalization explanation. 

After the completion of this work, we found other pa- 
pers where extended states have been seen in other ape- 
riodic lattices H]. However, the propagating modes in 
the FK model are somewhat unique as the aperiodicity of 
the FK model comes from a mismatch of the continuous 
sinusoidal potential with the periodicity of the lattice. 
The existence of propagating modes in the pinned FK 
model suggests the possibility of extended states in other 
pinned systems such as pinned flux and vortex lattices. 
Finally, the FK model describes a variety of systems such 
as adsorbed monolayers on substrates and charge density 
wave conductors [jl2|,p|, and we hope that our study will 
provide a new direction of research in complex aperiodic 
systems. 
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FIG. 1. Schematic diagram showing the equilibrium ball 
configurations in the incommensurate FK model for K large 
enough so that the effects of the order 0(K~( p+1 ^) are neg- 
ligible, (a) p = 2 and the resulting incommensurate phonon 
potential has two distinct levels which leads to the symbolic 
representation ...abaabaaaabaaa... of the lattice, (b) p = 3 
and there are three levels which generate the symbol sequence 
...acaacabbacabb.... By decimating the sites c out, one obtains 
the renormalized lattice ...aaaabbaabb... with dimer (aa) and 
tetramer (aaaa) "defects". 

FIG. 2. Phase diagram in the K — A parameter plane for 
E = E max and a = (3 — V5)/2. Here r = 2 arctan( K )/n 
and t = 2arctan(A)/7r. The dashed line corresponds to the 
standard FK limit. The smooth curves show the parameter 
values for the existence of stepon modes with k = fl = 0. 
The shaded part represents the phase where the eigenfunc- 
tions are exponentially localized. In the light shaded regime 
(K < K c ), the fluctuations are self-similar and fall into the 
universality class of the Harper equation whereas in the dark 
shaded regime (K > K c ) the fluctuations appear irregular. 
The unshaded part for K < K c corresponds to smooth ex- 
tended eigenfunctions while for K > K c the unshaded regime 
describes the critical phase. 

FIG. 3. The stepon solution at the first degeneracy point 
r = t = 0.67052674.... Here 6 = {na-a/2}. The stepon mode 
has its discontinuities at the same locations as the underlying 
phonon potential. 
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Ketoja and Satija, Fig. 1 




Fig3 Ketoja and Satija 
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